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1 Abstract 

This contribution deals with the creation of numerical models for the simulation of 
the dynamic characteristics of fractional-order control systems and their comparison 
with analytical models. We give the results of the comparison of dynamic properties 
in fractional- and integer-order systems with a controller, designed for an integer- 
order system as the best approximation to given fractional-order system. Other 
open questions are pointed out, which should be answered in this area of research. 

2 Introduction 

The standard control systems used so far were all considered as integer-order sys- 
tems, regardless of the reality. In their analysis and design, the Laplace transform 
was used heavily for simplicity. Because of the higher complexity and the absence of 
adequate mathematical tools, fractional-order dynamical systems were only treated 
marginally in the theory and practice of control systems, e.g. ||, |j. Their 
analysis requires familiarity of work with fractional-order derivatives and integrals 

ii, II, mum- 

By removing the restrictions to integer-order systems it is possible to obtain 
systems whose properties are a combination of systems of the closest integer-order, 
but also intermediate types of systems, which broadens the class of the systems 
considerably fl]]. With different fractional-order systems the notions arise such as 
weak or strong integrator or differentiator, weak or strong fractional-type pole, or 
zero, with interesting contribution to the dynamics of the system (stability, phase 
shift etc.), as some properties are retained, others are eliminated. A fractional-order 
system combines some characteristics of systems of the order N and (N + 1). Hence 
by changing the order as a real and not only integer value we have more possibilities 
for an adjustment of the roots of the characteristic equation according to special 
requirements. 

In this contribution we will analyze dynamic properties of systems in the time 
domain with an emphasis on the numerical methods of simulation of fractional-order 
systems. We will point out problems of inadequate approximation of fractional-order 
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systems with integer-order systems and the differences in dynamic properties of such 
systems in closed control systems with controller. 

3 Definition of the fractional-order control system 

For the definition of the control system we consider a simple unity feedback control 
system shown in Fig.l. G s (s) denotes the transfer function of the control system 
which is either integer-type (Gl(s)) or more generally fractional-type (G{(s)) and 
G r (s) is the transfer function of the controller, also either integer-type (G l r (s)) or 
fractional-type (G/(s)). 



W{s) + ^E(s) 


G r (s) 


U(s) 


G s (s) 


Y(s) 























Figure 1: Simple unity feedback control system 

For later purposes consider a fractional-order controlled system, which represents 
our real system, with the transfer function 



1 



1) 



a 2 s a + a\sP + a 

where a and (5 are in general real (a > (3). In the simulation, the coefficient values 



a2 = 0.8, a\ = 0.5, ao = 1, a = 2.2, (3 = 0.9 were chosen ]10|, |12j. To the fractional- 
type transfer function ([l]) there corresponds, in time domain, the fractional-order 
differential equation 



(2) 



a 2 y( a \t) + a iy W(t) + a y(t) = u(t) 
with initial conditions y^(0) = and y(0) = 0. 

4 Numerical and analytical computation of the 
unit step response of a fractional-order system 

For the numerical calculation of the unit step response of the fractional-order system 
(||) we employ, for the approximation of the fractional derivatives in equation (||), 
the relation (|3|) with "short memory" principle, formulated in [p~l] ] 



N(t) 

y (a) (t)^ {t - L) D?y{t) = h- a ^b jy (t-jh), 

j=0 



(3) 



where L is "memory length", h is time step, 

N(t) = min 

[z] is the integer part z, 
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^• = (-l) J (jJ (4) 

where (") is binomial coefficient. To calculate bj it is convenient to use the following 
recurrent relation 

6 = 1, bj = (1-1±^)6 J -_ 1 (5) 

It follows from the estimates derived in ][ll]] that in our case the normed error of 
such approximation is 



whence we have the following constraint for the choice of "memory length" L: 

L > 2 1 (7) 
% r («) 

where 5q is the maximum admissible normalized error and T(a) is the Gamma 
function. 

By using the relation (|3]) we can approximate the differential equation (^) a 
different mode. Our approximation [10, 12] of equation (|^) in discrete time steps 
t m (m = 2,3, ...) has the following form 

m m 

a 2 h~ a bjVm-j + aih^ 13 ^ c jVm-j + a y m = u m (8) 

j=0 j=0 

or 

m m 

a 2 h~ a (b y m + bjVm-j) + ai/i _/3 (c y m + Y c iVm- 3 ) + a y m = u m (9) 
j=i 3=1 

From the approximation (|9|) we can derive [JTC , 12], the following explicit recurrent 
relation for the calculation of the values y m (m = 2,3, ...) 

m m 

u m ~ a 2 h~ a bjy m -j - a x hr^ £ Cjy m -j 

a 2 h a bQ + aih Pcq + a 

with yo = 0, yi = 0, uq = and u m = 1 for m = 1, 2, .... 

This algorithm does not require iterational calculations, in contrast to the pro- 
cedure given in Q. 

For the analytical calculation of the unit step response of fractional-order systems 
(||) we apply the analytical form of the impulse response of such system ||. By 
integrating the impulse response of such system we obtain the following analytical 
form of the unit step response of such a fractional-order system 

y(t) - — V (~^ m (^] m M m+ V E (m) (-^L t a - f3 ) (11) 

where E\„(z) is the Mittag-Leffler function in two parameters, 

sffiM'I^M-gsnwnFbr- (fe = ' 1 ' 2 "" ) (12) 
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For the calculation of the Gamma function we used the algorithm from [jTj. 
The calculation of the unit step response directly according to (|IT| ) is numerically 
difficult. For the system (Q), at least 20 terms of the series were needed along with 
the extended precision of the real variables. To employ (11) further it is necessary 
to consider asymptotic methods of calculations. 




Fig. 2 Comparisom of unit sptep responses (10), (11) 

Fig. 2 shows a graphical comparison of numerically and analytically computed 
unit step responses of the system (||) with time step 0.1 seconds. The calculations 
show a good agreement of both methods. The differences depend on the size of the 
time step and on the method of approximation of the value y\. 

The value y± = in ( |io|) was calculated using the initial conditions for differential 
equation (|2|) and relation (||) with m = 1 



yi 



08) 



yoh h 13 



b h-1 3 



(13) 



Better results (see Fig. 3) we can obtain by approximation the value y\ from 
initial condition yo directly with equation (|l0| ) for m = 1. The advantage of such 
approximation is more evident for differential equations (|l7|), ([lS|) and ( |23| ) with 
derivatives of the function u(t) (unit step). 




L.20 2 .tO t .SO 



Fig.3 Comparisom of unit sptep responses (10), (11) 
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5 Approximation of fractional-order system 
with integer-order system 



In integer- type linear systems, higher-order systems usually are approximated - un- 
der certain conditions || - with second-order systems for simplicity. So far, however, 
fractional-order systems usually have not been considered. Regardless of the real- 
ity, their approximation is performed with close integer-type systems, albeit with 
unfavorable consequences. 

Let us approximate the fractional- type system (||), with the coefficients as given 
there, with the integer-type system of second order 

4y"(t) + a\y'(t) + ajy(t) = u(t) (14) 

The coefficients a\ of the integer-type system take on the values a l 2 = 
0.7414, a\ = 0.2313, a = 1 under the condition 

rain (15) 



m * 2 

£(*/-*!■) 

7=0 




The comparison of the unit step responses of both systems is shown in Fig. 4. It 
should be noted that a generally better approximation was achieved for a,/3<2. 

6 Design of the regulator 

For the system (O) we choose as an approximation to the system (|2|) the integer- 
order regulator [|l0| , [l~2|| with the transfer function 

G i r (s) = K + T d s (16) 

The regulator will be designed so that a unit step at the input of the closed 
regulation system in Fig.l will induce at the output an oscillatory unit step response 
with stability measure St = 2 and damping measure Tl = 0.4. Then the coefficients 
for ( |l~6| ) take on the values K = 20.5 and = 2.7343 . The differential equation of 
the integer-order closed regulation system has the form 

4y"(t) + (a\ + T d )y'(t) + (a + K)y(t) = Kw(t) + T d w'(t) (17) 
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It follows from an analysis of the roots of its characteristic equation that the re- 
quirements for stability measure and damping measure are satisfied. The permanent 
regulation deviation is 4.6% . For the closed regulation system ( |i~7| ) it is very simple 
to derive the analytical form of the unit step response. In a manner mentioned above 
it is also possible to obtain an approximating recurrent relation for the numerical 
calculations 

m 

Kw m + Tdh' 1 £ d,-w m -j + a 2 hr 2 (2y m _ 1 - y m - 2 ) + («i + T d )hr x y m _ x 
Vm = a 2 /i" 2 + (ai + T d )h~i + (a + K) 

form = 1,2,... , y = 0, y_i=0, Wq = Q and w m = 1 for m = 1,2, ... 

The unit step response is plotted in Fig. 5. From the results of numerical sim- 
ulation it also follows that the regulation area within 10 seconds was 0.71 and the 
permanent regulation deviation is up to 5%. 
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Fig. 5 The unit step response (18) of the system (17) 



We now apply the integer-order regulator, designed for the approximated system 
(model of the real system), to the original fractional- type controlled system @. The 
differential equation of the closed fractional-order regulation system has the form 

a 2 y {a) (t) + T d y'{t) + aiyW it) + (a + K)y{t) = Kw(t) + T d w'(t) (19) 

However, the analytical analysis of this equation is not as simple as with the 
equation (|l~7|). We will calculate the unit step response numerically with the aid of 
an approximating recurrent relation derived as mentioned above 

m mm 

Kw m + Tjh- 1 £ djWm-j - a 2 h~ a £ bjy m -j - a x h~P £ CjVm-j + T d hr x y m _x 
= j^o j=i j=l 

Vm a 2 h~ a b + ai/i-^co + Tah- 1 + (a + K) 

(20) 

form = 1,2,... , yo = 0, wq = and w m = 1 for m = l,2, ... 

To compare the numerical calculations we will also derive an analytical form of 
the unit step response of the closed regulating system with a fractional-order con- 
trolled system and an integer-type regulator. We will obtain the form by integrating 
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the impulse characteristic from pi for the equation of type (|l9|), whence 
00 (-l) m fa + K^ m m ' - x fc 



y(t) = E 



m=0 



ml 



( 



£a(m+l)-/3ft 




b 2 



(21) 



with ^2 = 7^ and c 2 = ^ ■ 

Fig. 6 shows a graphical comparison of analytically (|2l|) and numerically (p 



computed unit step responses of the system ( |19[) with time step 0.1 seconds. We 
can see the better results of the numerical calculations by approximation the value 
yi directly with equation (p0|) (a) than by approximation with equation (|3|) (b). 




Fig.6 Comparisom of unit sptep responses (20a), (20b) and (21) 



A comparison of the unit step responses of the closed regulation system with 
integer-order system (|l7|), with the fractional-order system ( |l9|) is shown in Fig. 7. 




£ . 67 3.33 



Fig. 7 Comparisom of unit sptep responses (17), (19) 



It can be seen that the dynamic properties of the closed regulation system with a 
fractional-order system and an integer-order regulator designed for an integer-order 
system which is an approximation of the fractional-order system are considerably 
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worse than with the integer-order system and the integer-order regulator. The reg- 
ulation area during 10 seconds is greater by 50%, the system stabilizes later and 
has larger surplus oscillations. The system (19) is much more sensitive to changes 
in parameters. For example fi~0| , 12|, at the change of T d to value 1 the system ( |l~9|) 
is just behind the border of stability (Fig. 8), and with another decrease of T d it is 
already unstable, whereas the system (|17|) is stable. 




Fig. 8 Comparisom of unit sptep responses (17), (19) with changed T d 



Hence disregarding the fractional order of the original system, its approximation 
by an integer-type system of second order and an application of a regulator designed 
for the approximating system to the original fractional-order system is not adequate 
for our case a > 2. 

For the fractional-order system we now use the fractional-order regulator (l^, 12 ] 
with the transfer function 

Gl{s) = K + T d s s (22) 

The differential equation of the closed regulation system with a fractional-order 
system and a fractional-order regulator has the form 

a 2 y (a) (t) + a lV W (t) + T d y^ (t) + (a + K)y(t) = Kw(t) + T d (t) (23) 

In this case the fractional derivative 5 of the unit step w s (t) is no longer equal 
to zero for t > 0. The numerical computation of the unit step response will be done 
through an approximating recurrent relation derived as mentioned above 

m m 

Kw m + T d h~ 5 Y djW m -j - a 2 h~ a Y bjy m -j- 
j=0 j=l 



a 2 h- a b + aih-Pco + T d h- 5 d + (a + K) 

m m 

aih~P Y Cjy m -j ~ T d h~ 5 Y djU 



izL_ : ZzL_ ( 2 4) 



for m = 1, 2, ... , yo = , u>o = and w m = 1, m = 1, 2, ... 

The analytical form of the unit step response of a closed regulation system with 
fractional-order controlled system and fractional-order regulator can be obtained 
p0| , |i~2H by dividing equation (|23| ) into two parts, according to the right-hand side. 
By taking the derivative of the order 5 of the impulse characteristic for the second 
part, by subsequent integer-type integration of the impulse characteristic || of both 
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parts, and their adding together, we obtain the following resulting analytical form 
of the unit step response for equation of the type 



m=0 



k=0 




a + K 



j.a(m+l)—6k 
(- : E {m) 



b 2 

C2 



(_fLL/*-/*| + 
a-f3,a+/3m-5k+l\ Q 1 )' 



(25) 



a-f3,a+/3rn-5(k+l) 



with b 2 = and c 2 = §?■ . 

Fig. 9 shows a graphical comparison of numerically and analytically computed 
unit step responses of the system (Eg). 




Fig. 9 Comparisom of unit sptep responses (24), (25) of the system (23) 



In Fig. 10, a comparison of the unit step responses of the closed regulation system 
with a fractional-order regulated system and the integer-order regulator ( |l~9|) and 
the fractional-order regulator (p3[), with only two, experimentally found, parameters 
Td = 3.7343 and 5 = 1.15, is given. 




Fig. 10 Comparisom of unit sptep responses of the systems (19) and (23) 
It follows from the simulations (Fig. 10) that the dynamic properties of the 
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system (19) with an integer-order regulator (Fig. 7) improved with the use of a 
fractional-order regulator (p3|). Therefore in the following it is necesary to deal 
with the methods of synthesizing the structure and parameters of fractional-order 
regulators. 



7 Conclusion 

In the contribution we gave relations for simulation of dynamic properties of 
fractional-order systems, via numerical as well as analytical methods. Also presented 
were the results of some simulations. The calculations showed a good agreement of 
the new numerical method and analytical method of simulation of the fractional- 
order control systems. 

It was shown that the fractional-order system (|2|) can be well approximated by 
the second-order system ([l4|) for a, (3 < 2. The integer-order regulator designed for 
the integer-order system (|14|) can in this case be applied with good results also to 
the fractional-order system. A worse approximation is achieved with a or (3 > 2. An 
application of the integer-order regulator, designed for the integer-order system as 
an approximation to such a fractional-order system is inadequate and with a change 
of system or regulator parameters can lead to system instability. 

In the future it is desirable to deal with the asymptotic methods in the analytical 
method of computation of the dynamic characteristics of fractional-order systems. 
Further, it is necessary to deal with the methods of identification of the fractional 
order and of the parameters of such systems. The most important task will be the 
elaboration of methods of synthesis of the structures and parameters of regulators 
for such types of systems. 
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